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Abstract 

We introduce a new velocity selection criterion for fronts propagating into 
unstable and metastable states. We restrict these fronts to large finite 
intervals in the comoving frame of reference and require that their centers 
be insensitive to the locations of the ends of the finite intervals, exhibiting 
thus effectively an approximate translation invariance. Only one mono- 
tonic front has this behavior and its velocity is the one that is physically 
selected. We present analytic results in the case of piecewise parabolic 
potentials, and numerical results in other cases. 
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1 The selection principle 



In many systems rendered suddenly unstable, propagating fronts appear. The de- 
termination of the speed of front propagation into an unstable state has attracted 
attention recently since it cannot be achieved by simply solving the differential 
equation in the comoving frame of reference on a one-dimensional infinite do- 
main. Indeed, there are many such solutions on such a domain, even though 
the propagating front in practice always relaxes to a unique shape and velocity. 
Selection principles have thus been formulated to determine the observable front, 
without having to solve directly the partial differential equation of motion for a 
range of initial conditions. These principles have involved concepts of linear and 
nonlinear marginal stability, of structural stability and of causality [1] , and all of 
them try to deal with the puzzle of the reduction of the multiple solutions to 
the single observed one. All these selection principles examine the wave from the 
viewpoint of the moving front, the corresponding wave equation being reduced 
then to an ordinary differential equation involving the speed v of propagation. 
This speed has a unique value if the front is invading a metastable state, but not 
so for the case of invasion into an unstable state. The latter case has consequently 
attracted most of the attention. 

The selection principles mentioned above have been based on the study of 
Fisher's dimensionless equation, ut —u xx +f(u), on the interval (—00,00) with 
/(0) — f(l) — 0, the states u — and u — 1 being unstable and stable fixed points 
respectively. Aronson and Weinberger [2] have shown that sufficiently localized 
initial conditions will evolve into an observable front propagating with speed v*, 
provided / > on (0,1). This selected speed v* is the lowest speed for which the 
above partial differential equation admits a monotonic front joining the stable 
state u — 1 to the unstable state u — 0, and it satisfies v* > 2y / /'(0). Any 
monotonic traveling wave solution u(£) of Fisher's equation, with £ = x — v t 
being the coordinate in the comoving frame of reference, is a propagating front 
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with speed v > and satisfies the "steady state" ordinary differential equation 

d 2 u du „, , 

^ + ^ + /M = o, (i) 

with u(— oo) = 1, -u(oo) = 0. Note that any speed t> > t> * will give a monotonic 
front, though not the observed one. The selection criteria seek to determine this 
lowest speed v*, without solving the initial value problem, selecting thus one 
among the multiple possible fronts invading the unstable state. If the selected 
speed is determined by the linear order terms alone we have the pulled case 



v* = 2a//' (0) (linear marginal stability). If linear analysis fails, we have the 



pushed case v* > 2y/'(0) (nonlinear marginal stability). In the case of fronts 
invading a metastable state though, only one front is possible on the interval 
(— oo, oo), and hence the selection is automatic [3] . 

We shall adopt here a unifying approach, showing that there exists a single 
selection principle at work for both the metastable and unstable cases, and for 
the regimes of linear and nonlinear marginal stability alike. This approach ex- 
ploits the fact that the observed front is translationally invariant in the comoving 
frame of reference, even on a large finite interval, in the sense that its location 
is effectively independent of the ends of the interval. We solve then the steady 
state equation of motion on a large finite interval for an arbitrary given speed 
v, subject to the appropriate boundary conditions, obtaining a certain solution. 
This solution, however, will have approximate translational symmetry, thus be- 
coming a physically observable front, only for a certain value of v. It is this 
value v* of v that is experimentally observed. Thus the selected front is the one 
that is ef fectively translationally invariant on a large finite interval, in the 
comoving frame of reference. 

Indeed, let us assume that u(£) is the solution of Eq. (1) subject to the 
boundary conditions u(Li) = 1 and u(L 2 ) = 0, with L\ <C L 2 . There is only one 
such solution for a given velocity v. We define the continuous potential V(u), 
where f(u) = —dV/du. We multiply Eq. (1) with du/d^ and integrate from L\ 
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to L 2 , obtaining then 

v(o)-v(i)-XM + XW (2) 

with = du/d^. If is going to be a physically observable front on this 
large, but finite, interval, it will have to be essentially translationally invariant. 
This means that du/d^ will be effectively zero in the regions close to the bound- 
aries, u having reached a fixed point there. Consequently w(Li) and w(L 2 ) will 
tend to zero, while w 2 (£)df; will be finite and practically independent of Li 
and L 2 , as L\ — > — oo and L 2 — > oo. Hence the speed v of Eq. (2) becomes 
independent of the endpoints of the interval, acquiring a unique value. In other 
words, only the front with that particular speed can be translationally invariant. 

The requirement that the front be independent of the ends of the finite interval 
selects therefore the speed 

v* = (V(0)-V(l))/ [ L2 w\0dd, (3) 

with Li — > — oo and L 2 — > oo, as the speed of the physically observed front. 
Note that no distinction has been made here between metastable and unstable 
states. Indeed, given any particular speed v, we can find a front interpolating 
between the stable and the unstable or metastable state, as we shall demonstrate 
below, provided the solution is found on a finite interval. As the boundaries go 
to infinity, the value of the speed is restricted to v* and the front becomes the 
one corresponding to the speed of Eq. (3). 

2 Analytic examples 

We shall demonstrate the proposed selection principle through analytic and nu- 
merical work, both for the unstable and metastable cases. We shall adopt for our 
analytic work the continuous potential 

V(u) = a(l-4w 2 )/8 if 0<m<1/2, 

= i/(l-2u)(3-2u)/8 if u>l/2, (4) 
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where v > and a = 1 or — 1, depending on whether the state u = is unstable 
or metastable, respectively. The corresponding f(u), 

f(u) = au if < u < 1/2, 

= -u) if u>l/2, (5) 

is piecewise linear and leads to exact analytic solutions. A similar piecewise 
linear model, but with £ in the interval (—00,00), has been used for a discrete 
reaction-diffusion equation[4]. Piecewise parabolic potentials, always on the in- 
terval (—00,00), have been used in nucleation and crystallization problems as 
well [5]. 

Eq. (1) will be solved subject to the conditions u(Li) = 1, u(L 2 ) = 0, where 
Li <C <C L 2 . We shall assume that the midpoint of the front occurs at the point 
£0, with u(£o) = 1/2, noting that both w(£) and du/d£, have to be continuous at 
that point. There are thus five boundary conditions that have to be satisfied, 
whereas the solution will involve five unknown parameters for any given value of 
v, namely £ an d two constants for each linear piece of f(u). We expect therefore 
a unique solution for each value of v. 

Indeed, the exact solution of Eq. (1) for the / of Eq. (5) is 

yq) 2 sinh[V^+4^(6 - Li)/2] J '-^-^ 



l r ^_ £V2 sinh[y^4^(£-L 2) /2] o<e<L2 , 
2 sinhfvV - 4a (£ - L 2 )/2] " " 



(6) 



where £ satisfies 

2v = Vv 2 + Av coth[vV + 4z/(£ - + Vv 2 - Aa coth[vV - 4a(£ - L 2 )/2}. 

(7) 

These relations hold even if v 2 < Aa. The solution of Eq. (7) gives £ as a 
function of the speed v. 

(a) Unstable case : a — 1. 

(i) We examine the case v > 2 first. If we require <C £ ^ -^2, then Eq. (7) 



reduces to yv 2 — A = \fv 2 + Au — 2v. This equation has a real solution, v = v*= 



(z/+l) / y/2u — 2, provided v > 3. In other words, for that particular value of v the 
midpoint £o can be anywhere in the interval and cannot be determined, rendering 
thus the front effectively translationally invariant on the finite domain [Li,L 2 ]. 
The value v* is therefore the selected speed if v > 3. In fact, the departures of 
v from v* are miniscule and fall exponentially as the ends of the interval go to 
infinity. 

Since v * > 2 when v > 3, we have a pushed case. If however v < 3, Eq. (7) 
yields 

4o~^i+2 , (8) 

V u + 4z/ 

and the front sticks to the left boundary. 

(ii) Consider the case v < 2 next. Eq. (6) shows then that on the interval 
[£o, L 2 ] the solution u(£) is proportional to sin[\/4 — v 2 (£— L 2 )/2] and has multiple 
roots. Our front though has to be monotonic, with no overshooting occurring. 
Thus must become zero only on the boundary L 2 . Therefore the point r = L 2 
—2ix/\/A — v 2 , where the above sine becomes zero, must lie outside the interval 
[£o, L 2 ], and hence £ > r. In that case Eq. (7) yields 

7T 2 r-2v + Vv 2 + 4z/coth[ v / v 2 + 4z/(^ - Li)/2]-i 

Co — -^2 /, + /, = tan 



V4 - V4 - v 2 L V4 - 

(9) 

This exact equation holds for v < 2 and all values of v, and determines the 
location £ of the midpoint of the monotonic front. Thus Eqs. (7) and (9) 
determine fully £ for a given v. In fact, if v < 3 and v is just below 2, Eq. (9) 
yields 



Co w £ 2 - 27r/V4 - t; 2 - (v/^TT - 2)- 1 , (10) 

or equivalently 



v » ^4 - 4tt 2 /[L 2 - Co + (2 - v/^TT)- 1 "] 2 . 

Thus f has a plateau as a function of £o at the value v — 2, if v < 3. 

Fig. (1) shows f as a function of the midpoint £o f° r the cases v — 9 and u — 2, 
for the potential of Eq. (4) with a = 1. For z/ = 9, the selected speed is t> * = 2.5 
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(pushed case). We see that for v > 2.5 the front is located close to Li, while 
it shifts abruptly to L 2 when the speed becomes less than v*. When v is equal 
to v*, the graph has a plateau, indicating that £ is pretty much undetermined, 
the solution being thus effectively translationally invariant. For the case v — 2, 
on the other hand, the front is located at L x when v > 2, as indicated by Eq. 
(8), while for v just below 2 the location of the front shifts abruptly to L 2 , as 
indicated by Eq. (10). We see thus the appearance of a plateau at v — 2, if 
v — 2. For that value of v the front's location is rather undetermined, indicating 
that the front has acquired effectively a translational invariance. The selected 
speed is thus v* — 2 (pulled case). 

Our analytic example indicates then that for v > 3 we have the pushed case, 
the selected speed being (v+ 1)/ \f1v — 2, while for v < 3 we have the pulled case, 
the selected speed being v* = 2. In both cases the selected speed corresponds to 
a plateau in the graph of v versus £o> due to the emergence of an approximate 
translational invariance of the front at that speed. 

(b) Metastable case : a = — 1. 

If v > 1 the stable state is at u — 1 and the metastable one at u — 0. The 
boundary conditions are once more u(L\) = 1, u(Li2) = 0, tt(£o) — 1/2, along 
with continuity of u and du/d^ at £o- For a given value of v there are then five 
boundary conditions and five unknown parameters, two constants for each linear 
piece plus £ - There is thus a single solution of Eq. (1) for any given value of 
v. This is the solution given by Eqs. (6) and (7), but with a — — 1. We can 
show that Eq. (7) reduces to v = v* = (v - l)/y/2u + 2 if L x < £o < ^2- In 
other words, for that particular value of v the midpoint £o can be anywhere in the 
interval, leaving the location of the front undetermined. The front acquires thus 
an effective translational invariance at that speed. This translational invariance 
is recognised in Fig. (2a), where v = 7, as a plateau in the graph of v versus £ 
at the speed v* = 1.5. 

Note that there is a solution to the metastable problem for any v on a finite 
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domain in the comoving frame of reference, but only for one v on the infinite 
domain. Indeed, if the metastable problem is solved on (—00, 00), then £ cannot 
be determined due to the exact translational symmetry, leaving us thus with five 
boundary conditions but only four unknown parameters, two for each linear piece. 
Consequently v will also have to be considered as a parameter to be determined, 
giving thus a solution only for a unique value of v[3], which is precisely the one 
selected by our selection principle. On a semi-infinite domain, on the other hand, 
a continuum of values is possible for v[6}. 

In the unstable case, exact translational symmetry on (—00, 00) reduces again 
the number of unknown parameters by one (£ ), but the boundary condition at 
u = is trivially satisfied due to the existence of two decaying exponentials, and 
thus we are left with four boundary conditions and four unknown parameters, 
two for each linear piece, for a given value of v. Solutions are thus possible for a 
continuum of values of v. 

3 Numerical examples 

We can demonstrate our selection principle numerically as well, for the case 

f(u)=u(b + u)(l-u)/b, (11) 

where for < b < 1 the states u = —b, u = and u = 1 are metastable, 
unstable and stable, respectively. In fact, it was this particular choice of f(u) 
that was used when the concepts of linear and nonlinear marginal stability were 
first proposed[7]. That study found that for 1 > b > 1/2 the selected speed 
for the front invading the unstable state is v* = 2, while for < b < 1/2 it is 
v* = (26+ \)j\p2b. For the front invading the metastable state the selected speed 
was found to be (1 — b)j\p2h. 

We have solved Eq. (1) numerically on a finite £ domain for the f(u) of Eq. 
(11), with 6 = 1/8, subject to the boundary conditions u(Li) = 1, u(L 2 ) = —b 
(metastable case). We found that the plot of v versus the characteristic point £ 
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of the front, where £o is defined through the relation «(£o) — 1/2, has a plateau 
at v — 1.75 (see Fig. 2b), indicating that at that speed the solution has become 
approximately translationally invariant on the finite £ domain. The speed on the 
plateau is precisely the one selected by the marginal stability criterion[7]. 

Furthermore, we have solved Eq. (1) numerically on a finite £ domain for the 
f(u) of Eq. (11), with b = 1 and b = 1/8, subject to the boundary conditions 
u(Li) = 1, u(L 2 ) = (unstable case). We find again that the plot of v versus 
£o, where tt(£o) = 1/2, has a plateau at v — 2 and v = 2.5, respectively (see Fig. 
3), indicating that the solution acquires effectively translational invariance there. 
These values are once again the ones known to be selected[7]. 

4 Concluding remarks 

We see then that requiring the solution to have approximate translational invari- 
ance on a finite interval in the comoving frame of reference results in the selection 
of a speed for the front in both the metastable and unstable cases. This speed 
is precisely the one given by marginal stability. A similar selection of a single 
velocity occurs when a cutoff is introduced, albeit on an infinite domain[8]. The 
importance of the translational invariance has also been noted in connection with 
the precursors of the propagating fronts [9]. Indeed, the selected solution is the 
only one with a legitimate translation mode in its stability spectrum. Thus one 
way of understanding the marginal stability of the selected solution is through 
requiring that the stability operator of a physically realizable solution possess a 
translation zero mode. 

We can adopt then a selection principle that reads "the selected front is the 
one that is approximately translationally invariant on a large finite interval, with 
respect to the comoving frame of reference" . This principle is very easy to imple- 
ment for both the metastable and unstable cases, especially numerically. Indeed, 
it suffices to solve the comoving frame equation on a large finite interval. For large 
speeds we expect the midpoint of the front to be close to the left boundary. As 
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the speed v is lowered, the midpoint suddenly moves to the right boundary. The 
speed v* at which this sudden move occurs is the speed selected by the physically 
observed front. 
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Figure Captions 



Figure 1: The speed v as a function of the midpoint £o of the front invading the 
unstable state, for the f(u) of Eq. (5), with L\ = —15 and L 2 = 15. The 
plateau is at v — 2.5 for v — 9, and at v — 2 for v — 2. All quantities are 
dimensionless. 

Figure 2: (a) The speed v as a function of the midpoint £ of the front invading the 
metastable state, for the f(u) of Eq. (5), with v — 7, L x — — 15 and 
L 2 = 15. The plateau is at v — 1.5. (b) The speed v as a function of the 
point £o of the front invading the metastable state, for the f(u) of Eq. (11), 
with b = 1/8, L\ = —15 and L 2 = 15. The plateau is at v — 1.75. 

Figure 3: The speed v as a function of the midpoint £o of the front invading the 
unstable state, for the f(u) of Eq. (11), with L\ = —15 and L 2 = 15. The 
plateau is at v = 2 for b = 1, and at v — 2.5 for b = 1/8. 
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